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The concept of duality has proved extremely powerful in extending our understand- 
ing in many areas of physics.^'^ Charge-vortex duality has been proposed^''' as a model 
to understand the superconductor to insulator transition (SIT) in disordered thin films^'* 
and Josephson junction arrays7'^'^ In this model, on the superconducting side, one has 
delocalized Cooper pairs but localized vortices; while on the insulating side, one has local- 
ized Cooper pairs but mobile vortices. Here we show a new experimental manifestation 
of this duality in the electron gas that forms at the interface between LaAlOa (LAO) and 
SrTiOa (STO).^"'^^'^^'^^'^'* The effect is due to the motion of vortices generated by the mag- 
netization dynamics of the ferromagnet that also forms at the same interface,^'*' ^^'^^ which 
results in an increase in resistance on the superconducting side of the transition, but an 

increase in conductance on the insulating side. 

The two dimensional SIT has been studied extensively^'*'^ as experimentally it is one of the few easily ac- 
cessible manifestations of a quantum phase transition. The two dimensional electron gas formed at the LAO-STO 
interface has been shown to undergo a SIT as the density of the carriers is tuned.'^' '''' In addition to superconduc- 
tivity, our earUer work showed evidence of ferromagnetism, which manifests itself as a hysteresis in measurements 

of the magnetoresistance (MR), Hall effect and superconducting phase boundary. Two other recent studies have 
confirmed the coexistence of superconductivity and magnetism in similar samples. The experimental evidence 
obtained so far'"*' indicates that the effect of the ferromagnetism on the electrical transport properties of the sys- 
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tern can be accounted for by considering only the magnetic field produced by the ferromagnet. As we show below, 
this close proximity of the ferromagnet and the superconductor allows for the observation of a novel demonstration 
of charge-vortex duality. 

The samples in this work had 10 unit cells (uc) of LAO grown by pulsed laser deposition on Ti02 terminated 
(001) STO single crystal substrates. Details of sample fabrication were discussed in an earlier publication^"* 
(see also Supplementary Information (SI)). The interaction between ferromagnetism and superconductivity can 
be seen in MR measurements in the superconducting state. Figures la and lb show the MR in a magnetic field 
parallel and perpendicular to the plane of the sample respectively, at T = 50 mK and Vg = 80 V, a gate voltage 
at which the sample is superconducting. Consider first the case of a parallel magnetic field, . As has been seen 
earlier for this sample in perpendicular magnetic fields, the parallel field MR is hysteretic due to the presence of 
ferromagnetic order (inset of Fig. la). Since is in the plane of the two-dimensional superconductor, the critical 
field is large {Hc\\ 1.2 T, see Fig. SI of the SI), so that one expects the MR to be small in this field range. In 
fact, the overall background MR is indeed quite small: most of the contribution to the MR is due to a set of twin 
peaks at ±12 mT. For clarity, we show in Fig. la the MR for only one sweep direction. As is increased from 
negative values, there is a small decrease in the resistance of the sample until H^^=0. The peak in the MR develops 
just as becomes positive for this field direction, but then rapidly dies out as iJy is increased father. The most 
surprising fact about the MR of Fig. la is that the amplitude of the peak depends on the sweep rate of iJy : the 
amplitude increases as the field sweep rate is increased. For very slow sweep rates, the peak is almost absent, and 
the amplitude of the overall MR is very small, as expected from our earlier discussion. 

The rate dependence of the MR can be understood if we consider the influence of the domain wall of the 
ferromagnet on the superconductor, as shown schematicaUy in Fig. 2a. In the simplest model, as H^^ is swept fi-om 
a large negative to a large positive value, the magnetization of the ferromagnetic layer reverses by forming a domain 
wall that propagates from one end of the sample to the other once i/y > 0. If the domain wall is a Bloch wall,^" the 
magnetization points out of the plane of the interface, and the conducting layer sees a magnetic field perpendicular 
to the interface that is localized near the domain wall, and travels along with it. We emphasize that for the real 
sample, the magnetization dynamics are likely far more comphcated, but the end result of a locahzed magnetic field 
perpendicular to the plane of the interface that moves as is swept is the same. At the field sweep rates in the 
experiment, the motion of the domains is quasi-static, a fact that can be confirmed by micromagnetic simulations. 
Consequently, the time dependence of the locahzed perpendicular magnetic field is determined only by the sweep 
rate of the external magnetic field. The moving perpendicular component of the field in turn generates moving 
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vortices in the superconductor. For our samples, which can be thought of as granular superconductors, moving 
vortices will result in a change in resistance^' as the vortices cross the weak Unks between superconducting grains 
(see section 6, Supplementary Information). This resistance change is proportional to the field sweep rate B and 
exponentially dependent on Ej and hence on the critical current 7c- One would therefore expect an increase in 
amplitude of the resistance peak with increasing field sweep rate at a specific gate voltage, and whose amplitude 
varies exponentially with at a specific sweep rate. This is exactly what is seen in our experiments (lower part of 
Fig. 4a, and Fig. 4c). 

Similar behaviour is also observed for the MR in a perpendicular external magnetic field H±_ (Fig. lb). Given 
the two-dimensional nature of the ferromagnet,^^ it is unlikely that the perpendicular magnetic fields appUed can 
cant the moments appreciably out of plane. However, a Bloch domain wall may form at the interface, and it is 
domains within this Bloch wall that are reoriented in an extemal magnetic field, as shown in Fig. 2b. Such behav- 
ior has been seen before in thin magnetic films.^^ (We stress again that this picture is highly simpUfied, and the 
real reversal process is Ukely much more compUcated.) H±_ itself leads to a large overall MR, with a minimum 
at zero field. The field from the ferromagnet gives rise to a sweep-rate dependent peak at H±^ ~ ± 15 mT, and 
a second, less prominent peak at H±^ ^ ± 48 mT, which is likely due to more complex magnetization dynamics 
in perpendicular field. Further evidence that the hysteresis and the dips are due to magnetization dynamics in the 
ferromagnet can be seen from the MR in crossed parallel and perpendicular fields. When one measures the back- 
ground MR in a perpendicular field in the presence of a large constant parallel field that aligns the magnetization 
of the ferromagnet, no hysteresis or rate dependence is observed (Fig. S2 of the SI). 

In the LAO/STO interface system, the superconducting transition can be tuned by VgP'^'^ Early studies of 
LAO/STO interfaces identified this as a SIT,'^ which has been studied extensively in the past in thin superconduct- 
ing films as a function of the film thickness or an applied magnetic field.^'* However, in the LAO/STO system, 
it appears difficult to tune the system deep into the insulating regime with a gate voltage.'^ Nevertheless, mea- 
surements of the current dependent differential resistance indicate that the transition is to a weakly insulating state 
where regions of superconductivity get increasingly isolated from one another as Vg is decreased (see Fig. S4 of 
the SI). Thus, as we noted above, it makes sense to model our sample as a granular superconductor consisting of 
a network of superconducting islands whose phase is coupled through the Josephson effect, parametrized by the 
Josephson coupUng energy Ej, and a Coulomb charging energy Ec, related to the energy cost of adding a Cooper 
pair to an island,.''-' (For a random network of islands, these quantities are averages over the network.) Changing 
Vg in our sample tunes the transition by modifying the ratio between Ej and E^. The superconducting and insulat- 
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ing regimes are duals of each other, and this duality is predicted to manifest itself through an interchange between 
specific measurable quantities. For example, the current / and voltage V are interchanged across the SIT, so 
that the I-V curve in the superconducting phase is similar in shape to the V-I curve in the insulating state. Such 
dual I-V characteristics have been observed in Josephson junction arrays and disordered thin fihns.^'^'^ 

In our samples, as noted above, it appears difficult to tune the system deep into the insulating regime, so that 
observing this signature of charge-vortex duality is difficult'^- (see the discussion in the SI). However, the inter- 
play between ferromagnetism and superconductivity at the LAO/STO interface results in a unique manifestation 
of charge- vortex duaUty associated with moving ferromagnetic domain walls that has not been observed in other 
systems. 

Figures Ic and Id demonstrate this new behaviour. Figure Ic shows the equivalent of the data of Fig. la but at a 
gate voltage of = -100 V, placing the sample on the insulating side of the transition. At the same magnetic fields 
at which a sweep-rate dependent peak was seen in Fig. la, a sweep-rate dependent dip is observed. As with the 
peak, the magnitude of the dip increases with increasing sweep rate. Similar behavior is observed in perpendicular 
magnetic field (Fig. Id). Note that in perpendicular field, one even observes two sets of dips (at ± 15 mT and ±48 
mT), mirroring the behavior seen in Fig. lb. In the charge-vortex duality model, current and voltage are switched 
in going across the SIT, so that conductance and resistance are also switched. Thus, where one observes peaks 
in resistance on the superconducting side, one should observe peaks in conductance (or dips in resistance) on the 
insulating side. This is exactly what is observed. As Vg is changed from -h80 V to -100 V, the peaks in the MR 
change to dips (Figs. 3a and b). A more striking graphical demonstration of the SIT can be seen if we plot the 
rate dependence of the resistance at the peaks or dips at different gate voltages normalized to their values at the 
highest sweep rate (Figs. 4a and b). We note that the resistance per square, Rs, at which the transition occurs is 
approximately 2.1 kO, less than the universal resistance value of /i/4e^ -^6.45 kfl expected for the SIT.^ 

The peaks in the MR on the superconducting side are due to the magnetic field associated with a moving 
nonplanar magnetization in the ferromagnet as the external magnetic field is swept. Given the fact that the position 
of the peaks and dips coincide for both H± and -ffy, it is clear that the dips also have the same origin. How 
would this field give rise to an increase in conductance on the insulating side? If the non-superconducting regime 
was simply a weakly localized metal, it would be hard to explain the dips in the MR that we observe (for a 
discussion see the SI). These dips can be explained if one models the insulating side of the transition as one where 
the Cooper pairs are localized on isolated superconducting islands, and a Cooper pair requires an energy Ec to 
transfer between two neighboring islands, even though Ec might be very small. Classically, the moving nonplanar 
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magnetization of the ferromagnet gives rise to a moving (i.e., time-dependent) field, which in turn gives rise to 
an electric field £h = —dA/dt through Faraday's law, , where A is the vector potential associated with the field. 
The voltage difference V = End generated between two neighboring superconducting islands separated by a 
distance d can exceed the charging energy E^., V > Ec/2e, so that the conductance of the system is momentarily 
increased. In contrast to the superconducting regime, this model predicts that the magnitude of the dip should 
depend exponentially on the sweep rate B (see Supplementary Information), a prediction that is borne out by the 
data (Fig. 4d). 

In sunnmary, the LAO/STO interface, with its unusual combination of superconductivity and ferromagnetism, 
provides a unique signature of charge-vortex duality in the superconductor-insulator transition. It would be in- 
teresting to see if similar signatures can be observed in artificially fabricated hybrid ferromagnet-superconductor 
structures. 
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Figure 1 : MR in the superconducting and insulating regimes, a, Parallel field MR as a function of different rates 
in the superconducting regime at = 80 V. Data for only one field sweep direction is shown for clarity. Arrow 
indicates the direction of field sweep. The inset shows the MR for forward and backward field sweep directions at 
a sweep rate of 300 f/T/s. b, Perpendicular field MR for fields swept from negative to positive values for different 
field sweep rates for Vg = 80 V. c, Parallel field MR for Vg = -100 V. d, Perpendicular field MR for -100 V. 
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Figure 2: Magnetization dynamics of the ferromagnet. a, Schematic of the system at different field values in 
parallel field. The top layer is the ferromagnet, the superconductor is shown through a gradient as extending some 
distance into the STO. The saturated configurations are shown at the two extremes when a high field is applied. 
Magnetization reversal by means of domain wall propagation in the ferromagnet for parallel field is depicted at a 
field where the peak in MR occurs. The domain wall induces vortices in the superconductor due to its perpendicular 
component of the field. It should be noted that the actual magnetization dynamics are more complicated in the real 
system, i.e. there are multiple domin walls moving in various directions, b, Schematic of the magnetization state 
of the system at different perpendicular fields. Due to the shape anisotropy of the system, the majority of the 
moments lie in plane, but the external magnetic field orients the direction of the perpendicular component of the 
magnetization of the domain wall. 
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Figure 3: SIT in the magnetoresistance. a, Parallel field MR for the fastest sweep rate (300 ^T/s) as a function 
of gate voltage, tuning the system from a superconducting to an insulating regime. The top and bottom panels are 
for insulating and superconducting regimes respectively. Charge-vortex duality manifests itself as the conversion 
of the peak in the superconducting regime to a dip in the insulating regime. The peak and dip occur at the external 
field value of ^ zL 12 mT. (There is an axis break on the y-axis in the bottom panel.) The maximum change in 
the resistance occurs for Vg = V, where the superconductivity is very weak, b, Similar behaviour in perpendicular 
field. Note that the fastest sweep rate in this case is 600 /iT/s. Additional structure in the perpendicular field at 
H± ^ ± 48 mT is due to the more complex magnetization dynamics in the perpendicular field. 
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Figure 4: Sweep rate dependence, a, The evolution of the resistance of the peaks and dips (normalized to the 
peak/dip value at the highest sweep rate) at H\\ ^ 12 mT with external field sweep rate as a function of the gate 
voltage for parallel field. The top and bottom panels are for insulating and superconducting regimes respectively. 
There is no saturation of the resistance of the peaks in the superconducting regime for the parallel field as the only 
contribution to pair destruction is from the field due to the domain wall, b, Similar plot for perpendicular field. 
The resistance of the peaks quickly saturates in the superconducting regime due to the pair breaking effect of the 
external field, unlike in the parallel field case. However, the resistance of the dips on the insulating side does not 
saturate, c, A plot of log{Rpeak — Rnopeak) vs the critical current, Ic, for the three fastest sweep rates in external 
field. Rnopeak is defined as the resistance at the same field value as Rpeak but on the field sweep curve in the 
opposite direction. The three lines are the fits to the data at the three different sweep rates, d, Rdip — Rnodip as 
a function of the parallel field sweep rate. The data for each gate voltage, Vg = -30, -40, -60, -80, and -100 V are 
fitted to an exponentially decaying function. 
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Supplementary Information for Evidence for charge- vortex 
duality at the LaAlOs/SrTiOs interface 

1 Experimental Techniques 

The samples in this work had 10 unit cells (uc) of LaAlOs (LAO) grown by pulsed laser deposition on Ti02 termi- 
nated (001) SrTiOa (STO) single crystal substrates.''^ The electrical measurements were performed on a Hall bar 
defined by photoUthography and etched using argon ion milling. Details of the film preparation, characterization 
and sample fabrication have been discussed in detail in prior publications.''^'^ The samples were measured in 
an Oxford dilution refrigerator with a base temperature of 15 mK. This refrigerator was equipped with a 2-axis 
magnet so that a magnetic field could be applied both perpendicular and parallel to the LAO/STO interface. A gate 
voltage Vg was applied to the back of the 500 /itm thick STO substrate. 

Transport measurements were made by a standard ac lock-in detection method. The samples were current 
biased with an excitation current of = 10 nA at a frequency of 11 .3 Hz. For = 20 V and higher, this excitation 
current is much less than the critical current, Ic (lac < 0.17c)- The voltage signals were first amplified using a 
low noise homemade instrumentation amphfier and then measured using a PAR 124 analog lock-in amplifier or an 
EG&G 7260 digital lock-in amplifier. 
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2 Determination of film thickness and parallel critical field 

From the continuous mapping of vs H±_ one can estimate the superconducting coherence length, ^, as discussed 
in Ref. 3. We obtained a value for the coherence length of ^ ~ 70 nm at 1^ =80 V. In order to determine the 
thickness of the superconducting layer, we performed a similar measurement of the dependence of on ify . The 
system was biased at the mid-point of the resistive transition {Rs = 288 O @ = 80 V), under a proportional- 
integral-differential (PID) feedback circuit. The output of the PID was used to drive the mixing chamber heater 
while the field was ramped. This ensured that we always stayed at T^, defined as the mid-point of the resistive 
transition. A plot of Tc vs measured in this way is shown in Fig SI. As expected for a 2D superconductor, the 
dependence of on a parallel applied field is quadratic except at low fields where hysteresis due to the ferromagnet 
suppresses Tc. From a fit to this measurement at high fields we extracted the parallel critical field, H^w ~ 1.2 T. 
For a thin film superconductor, Hc\\ and the film thickness are related by the formula'* 

where $o = /i/2e is the superconducting flux quantum. 

From a knowledge of ^ and H^w , we determined the thickness of the superconducting film to be d = 13.6 nm, 
which is less than ^. 

3 Magnetoresistance in crossed parallel and perpendicular fields 

Since our system shows hysteresis due to the magnetization dynamics of the ferromagnet, it becomes inherently 
more complex to study the effect of an externally applied field to the system. In order to eliminate the hysteresis 
we applied a persistent parallel field, ify, which saturated the magnetization of the ferromagnet in one direction, 
and then performed perpendicular field magnetoresistance measurements. In Fig. S2 we show data in which 
perpendicular magnetoresistance is measured for various values of -ffy for Vg = 100 V at T = 50 mK. It can 
be seen that for Hu =20 mT and higher, the sharp resistance peak in the magnetoresistance is absent and the 
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Figure SI: Parallel critical field Tc vs for Vg = 80 V. A parabolic fit (black curve) at high fields gives Hc\\ 
1.2 T. The thickness of the superconducting layer estimated from this value is d = 13.6 nm. 

magnetization of the ferromagnet is completely saturated. For = mT, the sharp peak is recovered. Note 
that the minimum of resistance and the field at which this minimum occurs both increase as the parallel field is 
increased. The additional pair breaking caused due to a constant parallel field gives rise to the small increase in 
resistance. The shifting of the zero of the magnetic field is due to a small misalignment of plane of the sample with 
respect to the axis of our split coil magnet. From the shift we calculate the misalignment to be 1.4°. 

In the insulating phase (for Vg = -40 V, say) the generation of flux due to domain wall motion is responsible 
for a rise in conductance. As is discussed in the main text and in detail below, a changing vector potential, dA/dt, 
generates an electric field, £h, which causes a potential, V, to develop between two adjacent islands. When this 
potential energy 2eV exceeds the charging energy of the islands, Ec, one gets conduction between the islands. In 
order to see that the dips in resistance that we see on the insulating side are indeed due to such a mechanism, we 
show R vs H± traces for Vg = -40 V, for = and 500 mT in Fig S3. Both the traces in Fig S3 are at the same 
sweep rate. By applying a parallel field of 500 mT, we have eliminated the effect of the magnetization dynamics in 
the ferromagnet. The sharp dips in resistance vanish completely, and only a small positive magnetoresistance due 
to weak localization remains. The shift in zero of the = 500 mT curve is due to the small misalignment of the 
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Figure S2: Crossed parallel and perpendicular field Rs vs H±^ for different values of parallel fields at Vg = 100 
V. For clarity field sweeps in only one direction are shown, with the arrow indicating the direction of field sweep. 

sample as discussed above. 

4 Insulating state I-V characteristics 

The I-V characteristics of the system were discussed in Ref. 3. The I-V curves in the insulating state were almost 
linear with no clear evidence of a Coulomb gap, indicating a weak insulating state. However, if dV/dl (Rs) is 
plotted as a function of an applied bias current Idc, a peak in dV/dl is observed at low bias. This is shown in Fig 
S4 for Vg = -30, -45, -60 and -100 V at T = 15 mK. For Vg = -30 and -45 V, a dip in dV/dl at zero bias signifies 
the presence of some superconductivity, however no evidence for that is seen in the temperature dependence of 
the resistance since dR/dT < at all temperatures in the temperature range of the current experiments. The dip 
in dV/ dl probably arises due to locaUzed superconducting islands weakly phase coupled by the Josephson effect. 
For Vg = -60 and -100 V, the system is in the insulating regime with a peak in dV/ dl at zero bias. 
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Figure S3: Charge-vortex duality Rs vs H± for — and 500 mT at Vg ~ -40 V. For both curves, the field is 
swept at the same rate. For — 500 mT, characteristic features due to the ferromagnet are absent. 

5 Weak localization effects in the insulating state 

It could be argued that the rise in resistance with decreasing temperature in the insulating regime is due to weak 
locaUzation effects in the electron gas. For a two-dimensional metal, the quantum interference contribution to the 
magnetoresistance including spin-orbit scattering can be written as^ 



ip is the digamma function, Hq = ^q/Stt^^, H2 = + (4/3)fl"so, = ^q/SttL^ and Hso = ^q/SttL^^. £ is 
the elastic mean free path, L^, is the electron phase coherence length and Lgo is the spin-orbit scattering length. 
Experimentally, the measured quantity is dR{H)/R = (AR/R){H) — (Ai?/i?)(0). Taking the appropriate 




(2) 



where 



hm = 7^(1/2 + HJH) - HHo/H), 



(3) 
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Figure S4: Insulating state I-V characteristics Rs vs I^c for different gate voltages at T = 15 mK. 



limits in Eqn. 1 above, one obtains: 
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where 



h{HO = ^{1/2 + HJH) - MH,/H). 



(5) 



It should be noted that Eqn. 3 does not depend on ifo> and consequently the elastic mean free path does not enter 
into the analysis. 

Figure S5 shows the 5R{H)/R for = -100 V at T = 50 and 400 mK. The curves at both the temperatures 
are reasonably well fit by Eq. 4. From the fitted curves one obtains at T = 50 mK, = 93 nm and Lso = 53 
nm, and at T = 400 mK, = 60 nm and Lgo = 42 nm. Using these values the WL contribution to the rise in 
resistance with decreasing temperature, SR{T, H — > 0), is estimated to be SR{T, iJ — > 0) ~ 200 f2 as T is reduced 
from 400 mK to 50 mK. However, the observed rise in resistance atVg = -100 V between these two temperatures 
is ^ 800 fl. Hence the increase in resistance cannot be solely accounted for by WL. 

The contribution to resistance due to electron-electron (e — e) interactions for a 2D metal is logarithmic in 



0.007 




- 



H|| = 100 mT - 



— fit 



-0.001 
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Figure S5: Weak localization analysis R{H) - R{0)/R{0) vs H± for Vg = -100 V at T 50 (red circles) and 
400 (cyan squares) mK in a constant parallel field, = 100 mT. The solid lines are the fits to the two curves 
using Eqn. 4. 



temperature and is given by^ 



where Tj-ef is some high temperature value for reference and a is a measure of the strength of the e — e interactions; 
typically a ^ 0(1). If we take Tref = 400 mK and calculate the e — e interaction contribution to the resistance at 
T = 50 mK, we get AR^^iT =50 mK) = ax 505 fl. Since the WL contribution - 200 n, and the total change in 
resistance ~ 800 fl, we can conclude that in our case a > 0(1). From the above analysis, it can be concluded that 
the major contribution to the rise in resistance with decreasing temperature in the insulating regime is due to e — e 
interactions. 




(6) 
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6 Thermal activation model for the resistance peaks/dips 

The model that we have for our sample is a random network of superconducting islands that are coupled with a 
characteristic Josephson energy Ej. The islands are small enough that there is a Coulomb penalty for transferring 
Cooper pairs from one island to another, characterized by the charging energy E^. Of course, for a random network 
of islands, there will be a distribution of Josephson couplings and charging energies, so Ej and Ec represent 
averages over the sample. 

6.1 Superconducting regime 

In this picture, on the superconducting side of the transition in the absence of motion of magnetic vortices, there 
is strong Josephson coupling between the islands {Ej » Ec), and the sample is in the zero resistance state. 
Any magnetic vortices present are localized in the interstitial regions between the islands. Motion of the magnetic 
vortices, which in our case is induced by the magnetization dynamics in the ferromagnet, will give rise to a finite 
resistance as the vortices cross the junctions between the superconducting islands. In order to understand the 
mechanism of the generation of this resistance, we use the phase shp model due to Langer and Ambegaokar (LA).^ 
There are corrections to this model, for example by McCumber and Halperin,' but the basic concepts are all that 
is essential for our analysis. 

To break down the problem, we consider two adjacent superconducting islands that are Josephson coupled to 
each other at two points, so that the interstitial region between them can enclose a magnetic vortex : the sample 
consists of a network of such interconnected islands. (More generally, one probably has multiple interconnected 
islands enclosing vortices, but this does not change the physical picture.) These two islands can be thought of as a 
dc SQUID, whose energy profile is given by 

U{<f>) = -E J cos <f> (7) 
where </> is the phase difference across the SQUID and Ej is the Josephson energy. Ej is proportional to the critical 
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current 7c of the SQUID. This defines a periodic potential where the minima of the potential correspond to integral 
values of the phase (/> = 2mr. At finite temperature, the system can be thermally activated from one potential 
minimum to an adjacent minimum which differs by a value Scj) = ±2^ over the energy barrier represented by Ej. 
Each such phase shp event will give rise to a voltage pulse according to the Josephson relation 

where the phase change is ±27r. In the absence of an external current, phase jumps in either direction are equally 
Ukely, hence the average voltage measured across the junction is still zero, although such phase shp events may 
give rise to voltage noise. 

In the presence of an external measuring current /, the SQUID's energy profile is modified to the "tilted 
washboard" potential 

U{cj>) = -Ejcoscj>-^Jcj>. (9) 

In this case, phase slip events corresponding to the system traveling "down" the washboard tilt are slightly more 
Ukely than those in the other direction, leading to a finite average voltage, and hence a finite resistance. However, 
the probabihty of such events is still quite small at low temperatures if / << Ic- 

LA's original paper focused on phase slips in a single weak Unk between two superconductors. They derived 
the rate at which phase slips occur in each direction 

j?± = ne-^^^/''-^ (10) 

where the attempt frequency O depends on microscopic parameters as well as external parameters of the circuit. 
Here AF^ is the free energy barrier for phase jumps in the the two directions^ 

AF^ ^2Ej±—. (11) 
2e 
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In the LA model, the application of a voltage V drives the generation of phase shps, and in steady state, the average 
voltage due to the generation of phase slips is equal to the applied voltage, giving the relation 

V = (27rfi/e)l^e-^^-'/'=«^ sinh(7rfiJ/2efcBT) (12) 

In the low current limit (/ <^ ekBT/h) the corresponding resistance is 

R = V/I= {■Kh^/e^kBT)ne-^^''''''^. (13) 

In our case, a phase slip is generated whenever a magnetic vortex Une crosses a weak Unk. The rate at which 
these phase slips are generated would then be proportional to the time dependence of the field generated by the 
magnetization of the ferromagnet. At the field sweep rates in our experiment, the magnetization change is essen- 
tially quasi-static, hence this time dependence is proportional to the extemal magnetic field sweep rate B, resulting 
in an additional factor in Eqn(13) proportional to B. The change in resistance due to the motion of vortices for the 
entire sample is an average of terms such as Eqn (13) with a distribution of Efs, but will be proportional to the 
magnetic field sweep rate B. The measured value of the critical current 7c of the sample will be determined by a 
parallel combination of random paths through the weak links connecting the superconducting islands. Neverthe- 
less, it is reasonable to assume that this measured critical current Ic is proportional to the mean Josephson energy 
Ej for the superconducting network. Thus, we can write the change in resistance due to the motion of magnetic 
vortices in the form 

Ai? = ABe-'^^'"/^"'^ (14) 

where Ic is the measured critical current, and A and a are numerical constants at a fixed temperature. Thus, if one 
plots In(Aii) as a function of Ic at a fixed temperature, one should obtain a straight line. Ic in this system can be 
varied by changing the gate voltage Vg. Figure 4(c) in the main text shows this plot. As can be seen, the expected 
exponential dependence is indeed observed for a nxmiber of different magnetic field sweep rates. The exponential 
factor a for the different sweep rates (which should be the same) matches to within about 20 %. 
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Equation (14) also predicts that the resistance change should be proportional to the magnetic field sweep rate 
B. If one looks at the parallel field data for Vg= and 20 V shown in Fig. 4(a) of the main text, one can see that 
the dependence on sweep rate is indeed approximately Unear. However, for Vg= 40, 60, 80 V, which are deeper in 
the superconducting regime, the dependence is clearly not Unear. In fact, deep in the superconducting regime, one 
might expect very little change in resistance, as the exponential factor involving 7c in Eqn (14) would suppress any 
resistance change. Thus the data at the three slowest sweep rates is not surprising. However, at higher sweep rates, 
the peak amplitude for these gate voltage values does increase. We do not know the reason for this, although one 
possibility is that at higher sweep rates, the effective Josephson coupling between islands is modified. We note that 
the data for the perpendicular field shown in Fig. 4(b) of the main text do not show a linear behavior, but in this 
case, the situation is comphcated by the fact that there is also a orbital contribution to the peak height due to the 
external perpendicular magnetic field. 

6.2 Insulating regime 

As discussed in the main text, on the insulating side of the transition, charge transport between islands is suppressed 
due to the charging energy Ec- If one considers thermal activation, the rate at which charge can be transferred is 
proportional to ^-^cI^bT ^ gnd as with the superconducting regime above, charge transfers in both directions are 
equally hkely, and there is no net average current. In the presence of a finite voltage difference ISV between the 
two islands, the rate is modified to p-i^c^'^e.i^VM^'BT ^ consequently favors charge transfer from the island 
at higher potential to the island at lower potential, resulting in a net average current between the islands. In our 
model, the potential results from the electric field that is generated by the moving magnetic field generated 
by the magnetization dynamics in the ferromagnet through Faraday's law. The time dependence of the changing 
magnetization, as we have pointed out above, is proportional to the sweep rate of the external magnetic field B, 
and hence ISV = 'yB, where 7 is a constant. Thus, the resistance of the dip should have a dependence on the 
magnetic sweep rate of the form 
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We note that this is a characteristically different dependence than that which is predicted and found in the super- 
conducting regime. 

Figure 4(c) in the main text shows a plot of the resistance change at the dip as a function of the sweep rate 
for different gate voltages on the insulating side: the lines are fits to the exponential form Eqn. (15). As can be 
seen, the fits are quite good. For Vg= -100, -80, -60 V, and -40 V the value of 7 is almost the same. For Vg= -30 
V, it is quite different. The time dependent magnetic field generates an electric field; the resulting voltage AV is 
proportional to the the electric field, but also to the spacing d between the islands. Hence 7 should also depend on 
d. As Vg is changed to bias the system closer to the transition from the insulating side, it is reasonable to expect 
that the effective spacing between islands decreases, resulting in a decrease in 7, as is observed in Fig. 4(c) of the 
main text. 

7 Alternative scenarios for the dip in the MR in the insulating regime 

As we cannot drive the system deep into the insulating regime by changing Vg, one might argue that the model 
of isolated superconducting islands is not apphcable, and one has instead a weakly localized metal. However, it 
is hard to explain the origin of the sweep-rate dependent dips in the MR for a metal. For example, if one had a 
metal, one possible scenario is that the moving magnetic field associated with the domain wall leads to momentary 
eddy current heating of the electron gas, resulting in a decrease in resistance as dRs /dT < in the insulating 
regime in our temperature range (Ref. 3). However, as the resistance of the film is large, the eddy current heating 
is expected to be small. In addition, a detailed comparison of the magnitude of the dips in the MR in comparison 
with the change in resistance as a function of temperature argues against this mechanism. Assuming that the 
moving magnetic field results in eddy current heating of the electron gas, one can estimate the increase in electron 
temperature by mapping the magnitude of the resistance change at a dip at the fastest sweep rate in parallel field, 
Rdip — -R(O) onto the -R vs T curve.^ Doing this, one obtains an estimated increase in electron temperature of 
AT 65 mK for Vg = -30 V, as shown in Fig. S6. However, similar heating effects should also be observed 
for gate voltages at which the sample is resistive, but shows a decrease in resistance , for e.g., Vg = -20 V, which 
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shows a peak in the MR (Fig. 3a of the main text). Indeed, eddy current heating in this Vg range should be larger, 
as the resistance is lower. For Vg — -20 V, Rpeak — ^(0) = 620 51; however, the i? vs T curve for Vg = -20 V 
does not show a corresponding increase in resistance (Fig. S6). Clearly for Vg = -20 V neither the sign nor the 
magnitude of AR is consistent with the eddy current heating scenario. We have attempted to come up with other 
potential mechanisms to explain the MR dips on the insulating regime, but there appears to be no other explanation 
that agrees with the experimental data other than the one that is discussed in the main text. 
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Figure S6: Eddy current heating i? vs T for Vg = -30 and -20 V. The vertical arrows indicate the size of AR 
{Rpeak/ dip ~ R{0)) Corresponding to the change in resistance at the peak or dip in the MR; the horizontal arrow 
shows the expected change in temperature if the dip in the MR is caused due to the eddy current heating. The 
explanation clearly does not work for the lower curve, as the change in resistance at the peak is larger than the 
change in resistance due to temperature. 
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